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Abstract 


We evaluate the Gutzwiller trace formula for the level density of classically 
chaotic systems by considering the level density in a bounded energy range 
and truncating its Fourier integral. This results in a limiting procedure which 
comprises a convergent semiclassical approximation to a well defined spectral 
quantity at each stage. We test this result on the spectrum of zeros of the 
Riemann zeta function, obtaining increasingly good approximations to the 
level density. The Fourier approach also explains the origin of the convergence 

problems encountered by the orbit truncation scheme. 
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The subject of “quantum chaos” comprises the study of quantum systems whose classical 
analogue is chaotic. One of the main tools is the Gutzwillcr trace formula, which provides a 
semiclassical approximation (SCA) to the fluctuating part of the quantum level density via 
an infinite sum over classical periodic orbits. A major problem with this sum is that it does 
not converge absolutely, making manipulations of the sum rather suspect, and moreover 
dependent on the ordering. 

An important question in this context is whether the Gutzwiller sum can be made to con¬ 
ditionally converge, by using an appropriate truncation procedure and possibly a reordering. 
Convergence in this sense means that one obtains better and better approximations to the 
level density, limited only by the SCA itself, as the truncation point goes to infinity. The 
most natural truncation scheme is to include only orbits whose period is bounded by some 
T max (we will call this scheme “orbit truncation”). Unfortunately, this has not had much 
numerical success. Gutzwiller in |1| used orbit truncation to evaluate the semiclassical level 
density for the anisotropic Kepler problem, and was able to reproduce the level density for 
low energies but not for higher ones — an incongruous result for an SCA. A mathematical 
model which is analogous to chaotic systems, the zeros of Riemann’s zeta function (see be¬ 
low), displays similar behavior: as one adds “periodic orbits” (prime numbers, in this case) 
to the sum, peaks in the level density become submerged in increasing oscillations j|]. In the 
hyperbola billiard there are indications that the sum as written does conditionally converge, 
however this is an exceptional system in which the Lyapunov exponent is larger than the 
Kolmogorov entropy 0. 

In this work we treat the problem of conditional convergence from a slightly different 
point of view. We consider d(e), the fluctuating part of the level density in a finite energy 
range, and define a purely quantum spectral quantity dq-{e) which converges to d{e) as 
T —> oo. dr(e) is constructed via a smooth truncation of the Fourier transform of d(e). 
We then use the Gutzwillcr sum to obtain the SCA to dr(s), which for any T involves an 
absolutely converging sum over periodic orbits. The limit T —» oo can now be taken safely, 
avoiding the convergence problems encountered by orbit truncation — the only remaining 
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errors are due to the SCA itself. This approach differs from previously suggested ones |4]-j6j 
in that it does not involve a resummation of the divergent tail of the sum or invoke the 
bootstrapping mechanism, but rather retains the intuitively appealing ordering by orbit 
period. 

We test this approach ,which we call “Fourier truncation”, on the spectrum of zeros of 
Riemann’s zeta function, and obtain very good numerical results. Moreover, we use the 
Fourier approach to analyze the orbit truncation scheme as well, and show that it leads to 
incorrect evaluation of the last few “frequency” components of the level density. Because of 
the exponential proliferation of orbits, these errors can become very large and overwhelm 
the resulting approximation. 

Our starting point is the fluctuating level density d(E) = ~ E m ) — d(E), where 

the {E m } are the set of eigenvalues of the Hamiltonian and d(E) is the smoothed level 
density (as given by e.g. Weyl’s form). We will confine ourselves to totally hyperbolic time- 
independent systems with a bounded phase space, in which case the E m are discrete and 
correspond to eigenenergies. 

We now select a finite energy range [E 0 ± AE/2] by means of a real and symmetric 
window function x(e), which satisfies the requirement x(|e| > AE/2) 0. A convenient 

form is the Gaussian, x( £ ) — exp(—e 2 /2 a 2 ). We then define the local fluctuating part of the 
level density d{e) = d(E 0 + e) x{ £ )i and expand it in a Fourier integral, 



(la) 


(lb) 


The value of d{r) using the true quantum level density is 
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where e m = E m — E 0 . Let us now define dr(e) as the approximation to d(e) obtained 


by truncating d{r) at r = T. We perform this truncation by means of another real and 
symmetric window function BV(t), which satisfies Wr(r) = 0 for |r| > T. Thus d(e) = 
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limr^oo dr(e), where 

dr{^) = 

From (H) dr{s) is given by 
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dr{e) = x(£m) w T {£m - e) , 
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where wr{£) = (27 rh) 1 f dr Wt(t) exp(— ire/h) is the inverse Fourier transform of Wt(t 
The simplest window one can take is the rectangular or “boxcar” window, 

R ( 1 for |t| < T, 

W|(r) = 

v 0 otherwise. v 11 ' 

Wj-(e) contains long oscillating tails, which fall off asymptotically as O^ 1 ). Since this 
falloff is only harmonic it may be possible for oscillations from many different delta peaks 
to interfere coherently to yield large spurious peaks. Thus it may be advisable to choose a 
different window function. 

The problem of estimating the spectrum of an infinite sequence from a finite subset is a 
common one in signal processing, and a variety of windows, with varying properties, have 
been developed 0 . One of the more popular is the Hanning window, 
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The largest sidelobe (i.e. tail oscillation) of w^(e) is ~ 0.025 the size of the main lobe, vs. 
~ 0.22 for the rectangular window. Moreover, the tails of w^(s) decay asymptotically as 
0(e~ 3 ), and so one need not worry about the cumulative effect of distant tails. 

In order to derive the SCA for dr{£) we start from the Gutzwiller sum, which approxi¬ 
mates the fluctuating part of the level density by a sum over classical periodic orbits ||, 


<*“(£) = l E t \AE) cos 
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Here the summation is over classical primitive periodic orbits p and repetitions r. S p , T p , and 
u p are the action, period, and Maslov index of the primitive orbit p, respectively. A p>r (E) = 
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T p /1 det(Mp —J)! 1 / 2 is the stability prefactor, where M p is the monodromy matrix for orbit p. 
We also restrict the energy range so that E 0 is large enough to be in the semiclassical domain, 
while A E is small enough so that A PiT (E) and T P (E) can be approximated by their values 
at E 0 , and the action can be replaced by its linearized form, S P (E) ~ S P (E 0 ) + eT p (E 0 ). In 
the following quantities without an energy argument refer to their values at E 0 . 

Next, we insert © into ([lT|) and exchange the order of the p summation and the e 
integration. The validity of this step constitutes our primary assumption, and is equivalent 
to ordering ([?;) by orbit period. We can then perform the inverse Fourier transform on each 
term of the sum to get 
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where y(r/h) is the Fourier transform of the window function y(e). For a Gaussian y (e) 


it is also a Gaussian, x(x) = exp[— (ox) 2 /2]. One can readily verify that (Q) is then 

absolutely convergent, since y effectively selects for each r only these orbits (primitive and 
otherwise) whose period is within a few h/a of r. 

Finally, inserting © into (|3j) and exchanging the summation and integration gives 
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where JF t (T, e) is a function which depends on the two windows we used, 
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Since Wt{t ) vanishes for |r| > T and y (r/h) falls off rapidly for \r/h\ > a -1 , d|?(e) only 
includes orbits whose period does not extend much beyond T + h/a. Thus for any finite T 
one is left with an essentially finite sum, which gives the SCA for the well defined quantum 
mechanical quantity dr{s)- In contrast, an orbit-truncated Gutzwiller sum is not in itself 
an SCA to anything. 
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Another way of defining the difference between the two truncation schemes is the fol¬ 
lowing. Suppose that the action is exactly linear in the energy (this is the case for e.g. 
the wavenumber level density for billiards). Then, approximating the whole spectrum using 
orbit-truncation corresponds to first taking a —> oo in ([]), and then T —> oo. However, 
from (jP) it is clear that one should first take the T —> oo limit to get the local spectrum, 
and only then extend to the whole real axis by taking the cr oo limit. Since the energy 
and time to infinity limits do not in general commute, the two truncation schemes can give 
differing results. 

We now test the Fourier truncation method for a particular model, the zeros of the 
Riemann zeta function on the critical line z — \ — iE. This is a mathematical, rather 
than a physical, model. However, it is well known that these zeros behave outwardly like 
the eigenvalues of a chaotic system without any symmetries (see e.g. @). Moreover, the 
fluctuating “level” density for the Riemann zeros on the critical line is exactly given by 

~ _ 00 lf)£f T) 

d ( E ) = ~ 2 EE 9 ,r /-r cos(rElogp) , (11) 

p r =1 Z7T VP 

where the p are the prime numbers. Eq. ( |TT| ) closely resembles the Gutzwiller trace formula 
(0), with h — 1, i/p — 0, T p — logp, S p = Elogp, and A p ^ r = — logp/27r v /p F (note that 
here the action is exactly linear in the energy). This “system” is often used in numerical 
investigations of the properties of various spectral sums, since it is a simple matter to 
generate millions of prime numbers, while Ending periodic orbits in real dynamical systems 
is a laborious process at best. 

Using (|TTD the equivalent expressions to (|D and H) for the Riemann zeros are 

d(r) = -EE r logP) ^ ol ° gP + *(r - r logp) e~ irE °^} (12a) 

dr{e) = —2 ffe EE^^rtrlog P,e) e '^‘-. (12b) 

We test these expressions numerically in Fig. [T]. We used 1000, 10,000, 100,000 and 5,761,456 
primes (i.e. all primes below 10 8 ), taking Eq = 69.5 (for historical reasons) and a Gaussian 
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x(e) with a = 2.5. T was chosen in each case so that all the Fourier elements considered 
could be evaluated to a good accuracy using the available primes (more about this later). 
The top 4 lines in Fig. [L] show the resulting approximate level density using a rectangular 


window VF-P(r). We clearly see the 5 peaks build up as T increases. As expected from (g), 


the amplitude of the intervening oscillations is not strongly dependent on T, but is more 


or less fixed by the form of Wr{r). This is also evident in the bottom 4 lines, where we 
show the approximate level density using (r). Here the oscillations are much smaller, at 


the price of an increase in the width of the peaks, as expected. We conclude that Fourier 
truncation appears to work very well for the Riemann zeros. 

It is instructive to compare Fig. |l] to the level density obtained by orbit truncation at the 
same place. In this case, as more primes are added to the sum peaks start to emerge, but 
are eventually obscured by increasing oscillations 0 . We can now understand this behavior 
by examining Fig. ^], where we show the Fourier spectrum of the truncated sum (|TT|) using 
the same primes as in Fig. |1|. It is apparent that, as more primes are added, a spectral peak 
builds up near r ~ T max = logp max (denoted by a dashed vertical line), where p max is the 
largest available prime. This peak eventually dominates the resulting level density. 

We can understand this peak in the following manner. For large positive r = T max we 
can approximate ( |12a|) by its second term, and also neglect all repetitions. We can then 
write dT ma x ( Tr na r) in the presence of orbit truncation in the form of another Fourier transform 
between the variables x and £, taken at £ = E 0 , 




d(T max , £) ® &*(£) 


(13) 


where O(x) is the step function, &(£) = |h(£) — i/27r£ is its Fourier transform, /u(x) is the 
weighted period density /u(x) = Yl p x e~^ x 5{x — logp), and ® denotes a convolution between 
quantities which are considered functions of £. Since &(£) contains long harmonic tails, 
^T ma x(^max) will also include ( 2ttE 0 )~ 1 times the DC (£ = 0) component of d(T max ). From 
the prime number theorem /i(x) ~ exp(x/2), and so for large enough T max we get 


max 
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Thus the discontinuous truncation “folds” into dr max (T max ) an exponentially increasing con¬ 
tribution from the DC component, giving rise to the peak in Fig. |2j. The inset of Fig. 2a 
shows d Tmax (T max ) vs. (O). and we see that the above prediction works quite well. 

Finally, we note that unlike in the Riemann case, typical dynamical systems also include 
Maslov phases. This changes the picture completely, as the DC term which is responsible for 
the divergent oscillations is then strongly suppressed. We therefore expect that the errors 
induced by the orbit truncation procedure would in such cases be less severe. 

In this work we proposed writing the SCA to the level density of a classically chaotic 
system by means of a limiting process, consisting at each stage of a convergent SCA to a well 
defined spectral quantity which itself converges to the level density in the limit. Numerical 
tests on the application of this method to the Riemann zeros yielded increasingly good 
approximations to the level density with no divergent oscillations. Moreover, we were able 
to understand the oscillations that appear in the conventional truncation approach as an 
edge effect resulting from the abrupt nature of the truncation. 

This work was funded by a grant from the MINERVA. 
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FIGURES 

FIG. 1. The level density computed using Fourier truncation, using (from bottom to top) 
1000, 10,000, 100,000 and 5,761,456 primes. Top 4 lines: Truncating using a rectangular (“boxcar”) 
window. Bottom 4 lines: Truncating using a Hanning window. Each line is shifted upwards by 2 
relative to the preceding one. The dashed lines indicate the location of the true Riemann zeros. 

FIG. 2. The Fourier spectrum of the orbit-truncated sum ([ll]), using 1000, 10,000, 100,000 
and 5,761,456 primes in subfigures a-d, respectively. The point r = T max is marked by a vertical 
dashed line. Inset — |^T max (T ma x)| vs. T max . The dashed line is the prediction of (|l4|). 
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